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Abstract 

Magnetic guiding fields in combination with energy dispersive semiconductor 
detectors have been employed already more than 50 years ago for in-beam in¬ 
ternal conversion electron spectroscopy. Even then it was recognized that effi¬ 
ciency modulations may appear as function of the electron energy, arising when 
electrons hit a baffle or miss the sensitive area of the detector. Current high pre¬ 
cision beta decay experiments of polarized neutrons with conceptional similar 
experimental devices resulted in a detailed study of the point spread function 
(PSF). The latter describes the radial probability distribution of mono-energetic 
electrons at the detector plane. Singularities occur as function of the radial de¬ 
tector coordinate which have been investigated and discussed by Sjue at al. 
(Rev. Scient. Instr. 86, 023102 (2015)), and Rubbers (arXiv:1501.05131vl 
[physics.ins-det]). In this comment a rather precise numerical representation 
of the PSF is presented and compared with approximations in the mentioned 
papers. 
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1. Introduction 

In two recent papers mm the radial spread of charged particles moving in 
a solenoidal magnetic guiding field has been investigated. The physical back- 
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ground behind this attempt is based on the fact that the distribution of the 
5 particles at a detector with finite radius is a potential source of systematic er¬ 
rors in high precision experiments. To the latter belongs the measurement of 
the beta asymmetry in the decay of polarized neutrons, see, e.g., Ref. [3]. In 
principle, this fact is known since long. In the late 60th of the last century 
nuclear spectroscopists developed solenoidal transport systems, equipped with 
10 Si(Li) detectors as energy dispersive elements, for in-beam internal conversion 
electron spectroscopy. It was pointed out already in one of the first publications 
in detail [4] that the phase relation between polar emission angle and the radial 
coordinate at a circular Si(Li) detector, or a baffle between target and detector, 
results in unwanted fluctuations of the transmission probability as function of 
15 the electron energy. Efficiency modulations on the 10 % level were reported, 
e.g., in Ref. [5]. In order to overcome this problem, the unwanted wiggles were 
averaged out by wobbling the magnetic field strength, see, e.g., also Ref. HE]. 

From the work of Sjue et al. [T] and Dubbers [2] it is now well known 
that such efficiency modulations originate from singularities which appear in 
20 the so-called mono-energetic point spread function (PSF) as function of the 
radius coordinate. These singularities have been mathematically treated in the 
mentioned papers by means of different approaches. The probability density at 
the detector plane is presented by Sjue et al. with the aid of an integral equation 
while explicitly by Dubbers, cf. Eq. (14) [T] and Eq. (12) [3], respectively. The 
25 ways to find solutions are quite different. Sjue et al. solve the integral equation 
numerically, while Dubbers presents and discusses analytical approximations. 

This contribution describes an alternative approach with which numerical 
solutions of arbitrary accuracy for the PSF can be obtained. It is based on the 
mathematically correct parameter representation of both, the radius coordinate 
30 R{cos9) and the probability density dP{cosd)/dR at the detector plane. Pa¬ 
rameter is cos 9 which intrinsically is a function of the polar emission angle 9 at 
the source. Notice that for rotational symmetry d{cos 9) is just proportional to 
the solid angle element did. In the next sectionfirst some mathematical details 
are described. Results will be compared in section with those presented in 
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Figure 1: (a) Impact loci curve at the detector plane for polar emission angle variations in the 
interval 1 > cos 9 > 0.2 with decrements Acos0 = 0.001. It is called in the text trajectory. 
Azimuthal emission angle (p = 0. (b) Scatter plot for cos 6 and p randomly distributed. 

Parameters of Ref. ^ taken: B = 0.326 Tesla, pc = 0.976 MeV, ro = p/{eB) = 0.01 m, zq = 
0.10 m, and zojrQ = 10. 


35 Ref. [U12] • The paper closes with conclusions in section 


2. Radial distribution at the detector plane 


For the sake of convenience the nomenclature of Sjue et ah, [T] will be 
adapted in the following. A right-handed coordinate system is defined with 
the magnetic field B coinciding with the z direction. The detector is placed 
40 in the (x,y) plane at a distance Zq from the origin of the coordinate system in 
which the point source is located. The charged particle starts with polar angle 0 
and azimuthal angle (p, the latter defined with respect to the y axis. The point 
of impact at the detector plane is given by DO Eq. (4)] 


r 


ro 



( 1 ) 


with ro = p/(c/ ■ B) the maximum projected orbital radius, p the momentum, 
45 and q the charge of the particle. Fig. (a) depicts the loci the radius vector 
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r/ro traverses according to Eq. Q if the cosine of the polar emission angle 
9 is decremented in steps Acosd = 0.001. This impact loci curve must be 
distinguished from the projected orbit of the electron on its way from the target 
to the detector and will be called in the following trajectory. 

50 The radius coordinate in the (x,y) plane reads ^ Eq. (12)] 

The differential probability dPc per normalized differential radius interval d^R/ro) 
is 


dPc{cos9) dPc{cosd) dcosd dPc/dcos9 

dR/ro dcosd d{R/ro) d{R/ro)/dcos9' 


For the example of an isotropically emitting source with dP{cos6)/dcos9 = 1 
and d{R/ro)/dcos6 calculated from Eq. ([^ one obtains after some algebraic 
55 manipulations 


dPc {cos 9) 
dR/ro 


rocos^d\/l — cos^d 

(1 - cos 2 d) zo cos + 2rocos3d sin ( 27 ^) 


(4) 


or normalized to the unit area 


dPc{cos9) 1 dPc(cosd) 

dA 2tt Rtq dR/ro 

It can be shown that Eq. ([^ agrees with Eq. (14) of Ref. [2]. 

Treating cos d as a free parameter, both the radial detector coordinate i?(cos d) 
and dPc{cos9)/dR can be calculated with Eq. ([^ and Q, respectively. If the 
60 parameter cosd is varied within the interval {1,0} one gets an impression how 
dPc/d{R/ro) evolves as function of R/ro- A corresponding parametric plot is 
depicted in Fig. (a). Shown are branches which start at R/ro = 0 and end 
again at R/ro = ^ for each closed trajectory at the detector plane depicted in 
Fig. a (a). However, what is wanted, is the sum of all individual contribution 
65 for a chosen normalized radius R/ro- This is the sum 

dP{R/ro) ( dPc(cosd|>) dPc {cos 9\<) \ 

^ \ dR/ro dR/ro j 
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Figure 2: (a) Parametric representation of the differential probability dPc/{dR/vQ) as function 
of the normalized radius coordinate R/tq. The parameter cos 6 has been varied in the interval 
{1, 0.27}, i.e. 73 % of the emitted intensity has been exhausted. For cos0 = 1 the curve starts 
at the origin, evolves with decreasing cos 0 to the first spike which corresponds to the largest 
radius of the inner trajectory in Fig. ^(a), returns to R/tq = 0 and evolves from there to the 
next spike, (b) Sum of 500 trajectories as function of the the normalized radius coordinate 
R/vq. Magnetic field and geometrical parameters the same as quoted in Fig. ^ 


with cos0|^ and cos0|^ the two solutions of the equation 



for the n*^ trajectory in the interval 


zq 

ro27r(n + Uf) 


< cos 9 < min 


■_2o_ 

-ro27r(n — 1 + nf) 



( 8 ) 


and Uf = floor(zo/('Co27r)). The lower limit of the summation m is the smallest 
integer for which > 1 holds with the maximum radius of the n*^ 

70 trajectory in the interval defined by the inequality Q. 


3. Results and Discussion 

In Figj^ (b) numerical results on the basis of Eq. are depicted. Totally 
500 trajectories at the detector plane of Fig. (a) have been taken into ac¬ 
count. The accuracy of the calculation at R/tq = 1.92 and R/rg = 0.6 has been 
75 estimated to be in the order of 1.8 %, and 0.6 %, respectively. The accuracy 
can be improved by extending the summation in Eq. ([^ over more trajectories. 
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Figure 3: (a) Random distribution as shown in Fig. ^ (b) analyzed at the detector plane 
with a radial resolution of 0.8 % or for the treated example of 80 /rm. (b) Comparison of the 
numerical true PSF of Dubbers [2] with the approximations presented in chapter 3 in red with 
the exact solution in black. 

However, it should be mentioned that the series expansion Eq. Q converges 
rather slowly. 

In the interval 1.92 < R/ro < 2 infinitely many spikes appear. It is not 
80 necessary to treat them in a mathematically exact manner. This fact is demon¬ 
strated in Fig. [^(a). A number of 10® impact points at the detector plane were 
generated by randomly distributing cos0 and ip in Eq. Q. A small sample 
is shown in Fig. (b). The generated distribution has been analyzed with an 
virtual detector of 0.8 % spatial resolution in radial direction, corresponding for 
85 the chosen example with rg = 1 cm to 80 /rm. It can clearly be seen that all 
the spikes in the mentioned interval result in a mean converging against 0.5. 
It should be mentioned that the same argument holds more or less for a finite 
beam spot size of the same order of magnitude. Fig. (a) is fully in accord 
with Fig. 4 of Sjue et al. [T]. Probably the procedure applied by the authors to 
90 generate their Fig. 4 was nothing else than what has been just described here. 

In Fig. (b) the exact results are compared with the approximations elabo¬ 
rated by Dubbers [5]. In due distance from the spikes and for outer trajectories 
in Fig. (a) with many revolutions, his approximation apparently seems to be 
rather good. For the important inner trajectory the approximation is rather 
95 poor, and even normalization is not preserved. 
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Finally the question should be addressed whether the considerations on the 
PSF of Ref. [HE] can be applied for field conhgurations other than homogeneous 
ones. It has been pointed out by Rubbers [H ch. 5] that in axially symmetric, 
continuously descending magnetic fields the general formulas remain valid after 
100 replacement of i? by i? • yjB{zQ)/B{{)) with i?(0) and B{zq) the magnetic helds 
on-axis at the source and detector position, respectively. It is certainly true 
that the underlying adiabatic invariance considerations makes some valuable 
statements on the movement of charged particles in inhomogeneous magnetic 
fields [SI p. 592f], see also Ref. [5|. However, whether in the fundamental 
105 Eq. (§ the phase Zq / (rp cos 9) of the cosine function can simply be replaced 
by a/B(0)R(zo) zo/[{Bp)p COS0], with {Bp)p = p/q, maintaining otherwise the 
functional dependence is questionable and needs to be discussed in terms of 
quantitative accuracy considerations. For more general magnetic field configu¬ 
rations it seems unlikely to find analytical solutions equivalent to Eq. <§ and 
110 Q and experimentalists would be well advised to investigate their instruments 
from the beginning by Monte-Carlo simulations performed with exact orbit cal¬ 
culations. 

4. Conclusions 

A method has been described with which mono-energetic point spread func- 
115 tions can be calculated with arbitrary accuracy for a homogeneous magnetic 
guiding field. It has been shown by Monte-Carlo simulations that a finite detec¬ 
tor resolution or a finite target spot size smear out the singularities for trajec¬ 
tories in the detector plane originating from polar emission angles approaching 
9 = -K12. The results of Sjue et al. [1] are fully in accord with the results 
120 obtained in this paper. Although Rubbers [5] presents the correct parameter 
representation for the probability density function, which is not explicitly quoted 
by Sjue et al. [1], his analytical approximations for the singularities appear for 
the innermost trajectories to be rather inaccurate with even a significant viola¬ 
tion of the normalization. 
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In any case, the subject addressed in Ref. mm is appealing and certainly 
beneficial for intuitional and educational purposes. 
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